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\  Abstract 

\J 

Numerical  solutions  of  the  Navier-Stokes  equa¬ 
tions  are  obtained  for  laminar  hypersonic  near¬ 
wake  flow  behind  cones  at  zero  angle  of  attack. 
Numerical  simulations  are  presented  for  a  cone 
with  a  7-degree  half-angle.  The  study  includes 
three  values  of  free-stream  Mach  number  (4 .32, 
6.32,  10.0)  and  four  values  of  Reynolds  num¬ 
ber  ranging  from  21500  to  172000  referenced  to 
base  diameter.  The  numerical  scheme  is  pat¬ 
terned  after  the  Alternating  Direction  Implicit 
method  of  Beam  and  Warming,  with  modifica¬ 
tion  to  include  finite-volume  discretization  and 
spatially  variate  time  step  for  accelerating  con¬ 
vergence  to  steady  state.  Comparisons  are  matte 
with  experimental  data  for  base  pressure,  length 
of  the  recirculating  region,  and  general  flow- field 
structure.  Effects  of  Mach  and  Reynolds  num¬ 
bers  on  the  laminar  near-wake  flowfield  struc¬ 
ture  are  reported. 

Introduction 

Considerable  amount  of  research,  both  theo¬ 
retical  and  experimental,  has  been  devoted  to¬ 
ward  understanding  the  complex  structure  of 
hypersonic  near-wake  flows.  Early  interest  fo¬ 
cused  on  predicting  base  pressure  or  base  drag. 
More  recently,  effort  has  been  directed  toward 
understanding  the  hypersonic  near-wake  flow 
structure  since  it  plays  an  important  role  in  the 
development  of  the  far  wake.  Application  of 
the  time-dependent,  compressible  full  Navier- 
Stokes  equations  to  the  hypersonic  near-wake 
flow  problem  is  very  attractive  because  it  has 
the  potential  for  complementing  the  simplified 


or  “component”  theoretical  models,  while  cir¬ 
cumventing  the  persistent  problems  of  support 
and  probe  interference  that  plague  the  exper¬ 
imental  work.  This  study  is  concerned  with 
the  effects  of  Mach  and  Reynolds  numbers  on 
the  structure  of  hypersonic  laminar  near-wake 
flows.  Laminar  flow’s  are  the  first  choice  be¬ 
cause  the  numerical  solutions  are  not  clouded 
by  the  uncertainties  of  turbulence  modeling. 
Furthermore,  the  present  study  establishes  a 
promising  approach  for  more  detailed  inves¬ 
tigation  of  the  near-wake  structure,  and  for 
evaluation  of  the  spurious  effects  introduced 
by  model  supports  and  probes  used  in  wind- 
tunnel  tests. 

The  brief  review  that  follows  concentrates 
on  laminar  wakes  created  by  cones  in  super¬ 
sonic  or  hypersonic  flight.  Crane1  has  reviewed 
theoretical  and  experimental  work  on  super¬ 
sonic  wakes,  and  he  traces  the  development  of 
theoretical  models  from  the  classical  paper  by 
Chapman2,  who  considered  a  simple  inviscid 
model  of  flow  about  a  “dead-air”  region  aft  of 
the  body  (see  Fig.  1).  In  reality,  viscous  effects 
are  important,  and  other  shocks  may  appear, 
such  as  the  lip  shock  associated  with  flow  re¬ 
compression  near  the  base  corner8,4.  Some  of 
the  flow  models  are  “component”  models  that 
break  down  the  wake  into  separate  but  inter¬ 
acting  regions  or  features,  each  more  tractable 
than  the  whole.  Thus,  Baum8  investigated  the 
boundary  layer  upstream  of  the  corner,  Ohren- 
berger  and  Baum5,  and  Weiss®  concentrated 
on  the  near  wake,  and  Zeiberg  and  Bleich7 
studied  the  region  beyond  the  neck  of  the  wake. 
Dash8  and  co-workers  have  extended  a  fam- 
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ily  of  component  models  to  wakes,  jets  and 
plumes. 

The  experimental  studies  of  near  wakes 
include  flight  tests,  ballistic-range  tests,  and 
wind-tunnel  tests.  The  flight  tests  reported 
by  Cassanto0  and  Bulmer10  include  base  pres¬ 
sure  and  heat  transfer  to  the  body.  Pho¬ 
tographs  taken  in  ballistic  ranges  by  Lyons 
et  a/11,  Koch12,  Levensteins  and  Krumins13 
and  others  provide  details  of  the  wake  ge¬ 
ometry  and  shock  location,  as  well  as  infor¬ 
mation  on  transition  to  turbulent  flow  (see 
Fig.  2).  Wind-tunnel  experiments  are  com¬ 
plicated  by  the  problem  of  support  interfer¬ 
ence.  Dayman14  pointed  out  that  even  the 
smallest  wire  supports  can  have  dramatic  ef¬ 
fects  on  the  structure  of  the  near  wake  This 
conclusion  is  supported  by  the  experience  of 
many  others  with  all  types  of  supports.  To 
circumvent  this  problem  a  number  of  experi¬ 
ments  have  been  performed  with  magnetically 
supported  models1 5  “  20 .  These  studies  are  at¬ 
tractive  from  the  standpoint  of  support  inter¬ 
ference,  but  they  are  not  free  from  tunnel  and 
probe  interference,  as  will  be  explained  when 
discussing  the  present  results,  which  are  com¬ 
pared  with  the  experiments  of  McLaughlin  et 
at 18  and  Blankson  and  Finston20. 

Correlations  of  various  features  of  the 
near  wake  have  been  supplied  by  several  au¬ 
thors  to  check  their  results  or  to  predict  needed 

I  A  A  |  An 

information,  such  as  base  pressure,  •  ’ 
wake  transition,23,24  the  end  of  the  recircula¬ 
tion  region25  or  the  location  of  the  separation 
line26'28. 

Numerical  simulation  by  means  of  Navier- 
Stokes  techniques  has  been  applied  to  projec¬ 
tile  wakes  by  Sahu  et  af20,  by  Deiwert30  and 
by  Thomas  et  a/31 .  Reklis  et  a/32  show  prelim¬ 
inary  results  for  the  wake  of  a  sharp  cone  in 
hypersonic  flight;  Reklis  and  Conti33  investi¬ 
gated  the  effects  of  Reynolds  number  at  Mach 
6.3.  The  present  investigation  extends  the  lat¬ 
ter  work  with  more  accurate  results,  and  a 
study  of  Mach  number  effects 


Mathematical  Formulation 

The  selected  ranges  of  Mach  and  Reynolds 
numbers  used  in  this  numerical  study  justify 
the  assumption  of  laminar  flow  throughout  the 
near  wake.  This  is  illustrated  in  Fig.  2,  where 
those  ranges  are  compared  with  various  tran¬ 
sition  data  and  correlations.  Hence  the  clas¬ 
sical  Navier-Stokes  equations  provide  a  realis¬ 
tic  mathematical  model  for  hypersonic  laminar 
near-wake  flow  problems. 

We  can  write  the  three-dimensional  un¬ 
steady  compressible  Navier-Stokes  equations 
in  strong  conservative  form  in  a  Cartesian  co¬ 
ordinate  system  as 

qt  fi  ^  gv  h*  =  1(Ti  +  ©y  +  fl2)  (la) 
where 

q=  \P>pu,pv,p'w,p(]T 

f  =  \pu,(p/')  +  pu2),puu,puw1  ii(p/-y  i- pt)]T 

g  =  \pv,pvu,(p/' y  t  pv2),pvw,v{ph  -  pc)f 

h-  \pw,pwu,pwv,  (p/')+  pe)]T. 

The  equation  of  state  for  a  thermally  and 
calorically  perfect  gas 

c2  =p/p  =  T  =  1(1-\){e-V*/2)  (16) 

is  required  to  form  a  closed  system.  In  equa¬ 
tions  (1)  u,v,w  represent  the  components  of 
the  velocity  V  along  the  coordinate  directions 
i ,y,z  respectively,  p  represents  the  pressure, 
p  the  density,  and  c  the  total  (internal  plus 
kinetic)  specific  energy.  Equation  (la)  is  ex¬ 
pressed  in  non-dimensional  form,  with  refer¬ 
ence  values  for  the  pressure,  density,  temper¬ 
ature  and  viscosity  taken  as  their  free-stream 
values;  the  reference  speed  is  the  free-stream 
speed  of  sound;  the  reference  energy  is  the 
square  of  the  free-stream  speed  of  sound;  the 
reference  length  is  the  diameter  of  the  base 
of  the  cone;  the  reference  time  is  constructed 
from  the  reference  speed  and  length,  and  the 
Reynolds  (Rt)  and  Mach  (M)  numbers  are 
based  on  these  reference  values. 
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To  enhance  the  numerical  accuracy  of 
boundary  conditions,  a  generalized  coordinate 
transformation  of  tire  governing  equations  is 
employed,  while  still  retaining  the  strong  con* 
servation  form.  The  boundary-conforming 
curvilinear  coordinate  system  is  subjected  to 
the  following  time-dependent  transformation, 

x,y,zj  <=>  £,r?,y,r  (2) 

£  £(*,!/, *.0 

n 1  =  n(x,y,z,/) 

?  =  f(x,y,z,t) 

T  —  t, 

yielding  the  following  system  of  equations: 

Qr  *  /.  *  Hr,  1  h  =  /?;'(A  ,  ©,  4  nf)  (3) 

where 

q=  Jq 

f  -  J  (£<q  +  fyg  +  £zh) 

{f=J  {mq  +  rizf  -+■  rjyg  +  rj,h) 

/»  =  J  (f«q  4  +  fyg  +  f,h) 

J  (S,t  +  +  s.n) 

0  =  J  (q*r  +  rjy©  +  q*fl) 
n  —  j  (?zr  -r  ^y©  +  ?*n), 

and  J  is  the  Jacobian  of  the  inverse  transfor¬ 
mation  (2) 

J  =  d(x,y,  z)/d(£, *?,?)• 

The  column  vectors  ,T,  ©,  17  can  be  conve¬ 
niently  decomposed  into  three  terms,  each  one 
containing  derivatives  of  the  flow  variables 
with  respect  to  only  one  coordinate  (indicated 
by  a  superscript  in  parentheses)  i.e., 

r -  f<{)  tn)  + 

©  =  ©<€)  +  ©(,,)  +  ©(0 

n-  ni()  -  n{r,]  +  i?(f). 


The  detailed  expansions  of  the  terms  in  the 
right-hand  sides  of  these  equations  are  given 
by  Thomas34. 

Equations  (3)  and  (lb)  are  the  full  gov¬ 
erning  equations;  for  the  present  application 
to  conical  vehicles  the  streamwise  coordinate  is 
£,  the  circumferential  coordinate  is  rj,  and  the 
cross-stream  coordinate  is  Since  the  flow  is 
axisymmetric,  derivatives  of  the  flow  variables 
with  respect  to  t]  are  zero.  However,  some  of 
the  metric  terms  give  non-zero  contributions 
in  this  direction,  giving  rise  to  a  source-like 
term  in  the  left-hand  side  of  (3).  This  term 
was  evaluated  analytically,  following  Nietubicz 
et  a/35.  Some  of  the  transport  terms  in  the 
right-hand  side  of  (3)  were  neglected  for  con¬ 
venience,  namely  ©,  J^^and  17^.  Equation 

(3)  then  reduces  to 

9r  -I-  /j  -r  /lf  +  G  = 

(4) 

where  G  represents  the  metric  source  term. 
In  the  “thin-layer”  approximation36  only  the 
term  on  the  right-hand  side  of  (4)  is  re¬ 
tained.  In  the  present  application  all  terms  of 

(4)  were  retained,  and  the  “cross-derivative” 
terms  /V’  and  17[f)  were  found  to  be  essential 
to  achieve  convergence.  This  is  attributed  to 
the  complex  character  of  the  flowfield,  which 
features  boundary  layers  running  in  the  (  and 
the  £  directions,  recirculating  flow,  and  mul¬ 
tiple  stagnation  points.  The  cross-derivative 
terms  are  cumbersome  to  handle  implicitly, 
and  were  instead  differenced  explicitly,  with¬ 
out  undue  ill  effects. 

Numerical  Method 

The  numerical  scheme  implemented  is  pat¬ 
terned  after  the  Alternate  Direction  Implicit 
(ADI)  method  of  Beam  and  Warming87.  Sev¬ 
eral  modifications  have  been  made  to  the  nu¬ 
merical  algorithm,  including  the  finite-volume 
discretization  of  Thomas  and  Lombard88  and 
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central-difference  formulas  that  are  second- 
order-accurate  spatially.  Conservative  artifi¬ 
cial  dissipation  terms  have  been  added  in  or¬ 
der  to  suppress  the  high-frequency  components 
that  appear  in  regions  of  severe  pressure  gradi¬ 
ent,  such  as  the  neighborhood  of  shock  waves. 
An  explicit  fourth-order  term  has  been  added 
to  the  right-hand  side  of  the  difference  equa¬ 
tions,  and  an  implicit  second-order  term  has 
been  added  to  the  left-hand  side  of  the  dis¬ 
cretized  equations.  Also,  to  accelerate  the  con¬ 
vergence  of  the  solution  to  steady  state,  a  spa¬ 
tially  variable  time  step  has  been  incorporated. 
In  this  procedure  each  cell  in  the  computa¬ 
tional  domain  is  advanced  to  steady  state  at 
the  maximum  local  time  step  allowed  by  the 
linear  stability  analysis.  On  grids  for  viscous 
flow,  where  grid  lines  are  highly  stretched  nor¬ 
mal  to  the  surface,  the  time  steps  can  vary 
by  several  orders  of  magnitude.  This  fact  dra¬ 
matically  affects  the  convergence  rate  of  the 
solution. 

Mesh  System 

A  body-fitted  coordinate  system  is  desirable, 
since  boundary  surfaces  in  physical  space  are 
mapped  onto  rectangular  surfaces  in  computa¬ 
tional  space,  to  allow  more  accurate  treatment 
of  boundary  conditions.  The  physical  space, 
which  is  L-  shaped,  requires  non-uniform  grid 
distribution  to  exploit  the  important  length 
scales  involved  in  this  flow  problem.  The  trans¬ 
formed  space  retains  the  same  topology  (L- 
shaped)  but  with  uniform  discretization. 

Boundary  Conditions 

The  boundary  conditions  are  as  follows  (re¬ 
fer  to  Fig.  3):  Adiabatic-wall  and  no-slip 
boundary  conditions  are  implemented  explic¬ 
itly  along  the  side  of  the  cone  (A-F  in  Fig.  3  ), 
and  on  its  base  (E-F).  Upstream  of  the  base, 
along  a  normal  to  the  axis,  ( A-B)  a  fixed  inflow 
condition  given  by  a  forebody  solution  of  the 
viscous  flow  over  a  sharp  cone,  obtained  in  the 
manner  described  by  Reklis  et  a/82.  Along  the 
axis  (E-D)  explicit  symmetry  boundary  condi¬ 


tions  are  used;  along  a  normal  to  the  axis,  at 
the  downstream  end  of  the  flowfield  (D-C)  an 
implicit  supersonic  outflow  condition;  and  on 
the  upper  side  of  the  flowfield  (B-C)  a  fixed 
shock  condition.  The  shock  is  simply  a  recti¬ 
linear  continuation  of  the  bow  shock  created  by 
the  sharp  cone,  and  is  unaffected  by  the  base 
flow  ahead  of  the  leading  edge  of  the  base  ex¬ 
pansion.  The  shock-expansion  interaction  and 
its  downstream  effect  on  the  flowfield  are  ne¬ 
glected  on  the  grounds  that  they  are  of  little 
consequence  to  the  near-wake  region. 

Results 

The  numerical  solutions  in  the  present  paper 
address  the  effects  of  Mach  and  Reynolds  num¬ 
bers  on  the  structure  of  laminar  near-wake 
flows  at  supersonic  and  hypersonic  speed.  The 
computations  simulate  the  laminar  near  wake 
behind  a  cone  having  a  half-angle  of  7  degrees, 
and  adiabatic  wall.  The  flow  conditions  inves¬ 
tigated  in  the  present  paper  cover  a  range  of 
Mach  and  Reynolds  numbers  as  given  in  Table 
1. 

Table  I 

Flow  Condition  Matrix 


Red 

M=4.32  jM=6  32  M=10.0 

21,500 

X 

43,000 

X  X 

86,000 

XXX 

172,000 

X  X 

The  numerical  solutions  at  Mach  num¬ 
bers  of  4.32  and  6.32  are  compared  with  the 
experimental  data  of  McLaughlin  et  all 8  and 
Blankson  and  Finston  20  respectively,  as  dis¬ 
cussed  below.  The  computational  grid  used  in 
all  present  calculations  is  a  92x58  mesh,  with 
92  grid  lines  in  the  axial  direction  and  58  grid 
lines  in  the  direction  normal  to  the  axis.  Non- 
uniform  grid  distribution  was  used  in  both  di¬ 
rections.  The  mesh  system  extends  about  5 


diameters  downstream  of  the  base.  A  typical 
mesh  system  used  in  the  present  calculations 
is  shown  on  Fig.  4. 

Related  Experiments 

Experimental  investigation  of  the  near  wake 
of  models  suspended  in  wind  tunnels  is  ham¬ 
pered  by  significant  interference  from  the  sus¬ 
pension  system,  as  indicated  by  Dayman14; 
on  the  other  hand,  the  information  provided 
by  free-flight  experiments,  although  free  from 
support  interference,  is  severely  limited  by 
the  experimental  environment.  In  an  ef¬ 
fort  to  circumvent  the  support  interference 
problem,  a  few  investigators  have  resorted 
to  magnetic  suspension  of  models,  notably 
Murman16,  McLaughlinel  a/18,  and  Blankson 
and  Finston21'.  Data  from  the  latter  two  in¬ 
vestigations,  taken  in  the  same  wind  tunnel 
at  different  Mach  numbers,  were  used  in  the 
present  study  to  validate  the  numerical  simula¬ 
tions.  These  data  consist  mainly  of  Pitot  probe 
measurements  and  recovery  temperature  mea¬ 
surements.  McLaughlin  also  measured  static 
pressure  with  two  types  of  probe.  All  these 
data  are  essentially  free  from  support  interfer¬ 
ence,  but  are  apparently  affected  by  extrane¬ 
ous  tunnel  effects.  For  instance,  the  Blank- 
son  data  were  taken  in  an  accelerating  stream 
(nominally  at  Mach  6.32)  that  ranged  from 
Mach  6.45  to  6.65  along  the  wake  under  study. 
The  angle  of  the  bow  shock,  deduced  from  the 
Pitot  data,  is  10.8  degrees  near  the  base  of  the 
model.  This  value  is  low  compared  with  the 
theoretical  value  of  11.6  degrees  for  inviscid 
conical  flow,  and  the  numerical  result  of  12.0 
degrees  obtained  in  the  present  study.  These 
effects  hamper  direct  comparison  of  presure 
and  temperature  profiles.  Therefore,  quantita¬ 
tive  comparisons  are  limited  to  overall  quanti¬ 
ties  such  as  base  pressure  and  location  of  the 
rear  stagnation  point.  Qualitative  comparison 
of  wake  feature s,  as  well  as  profiles  of  Pitot 
pressure  and  total  temperature,  which  were 
measured  experimentally,  show  good  agree¬ 
ment  with  the  present  numerical  simulations. 


Accuracy 

The  choice  of  a  proper  grid  that  resolves  all 
the  essential  physical  features  of  a  given  flow 
problem  is  an  important  issue  that  needs  to 
be  addressed  in  any  numerical  simulation.  A 
flowfield  simulation  with  an  excessively  large 
number  of  mesh  points  requires  unnecessary 
computing  time,  whereas  an  insufficient  num¬ 
ber  of  mesh  points  may  degrade  the  accuracy 
significantly.  The  effect  of  grid  resolution  on 
the  numerical  accuracy  has  been  studied  in  the 
present  work.  The  Blankson20  experimental 
data  at  Mach  number  of  6.32  and  Reynolds 
number  of  86000  has  been  selected  as  a  bench 
mark.  In  the  computations,  58  mesh  points 
were  distributed  non  uniformly  in  the  radial 
direction,  and  this  number  was  held  fixed.  The 
axial  mesh  was  refined  successively  to  yield 
three  mesh  systems  of  46  x.  58,  92  x  58,  and 
135  x  58  points.  It  was  found  that  the  location 
of  the  rear  stagnation  point  is  sensitive  to  low 
streamwise  mesh  resolution. 

Table  2 

Effect,  of  Mesh  Resolution  on  Location 
of  Rear  Stagnation  Point 


Table  2  shows  the  axial  location  of  the 
rear  stagnation  point  for  the  various  grids.  It 
is  clear  that  the  46  x  58  mesh  is  not  suffi¬ 
cient  to  resolve  the  near  wake  structure  and 
leads  to  significant  error.  The  solution  for  the 
92  x  58  mesh  compares  well  with  the  experi¬ 
mental  data  of  Blankson  20  (see  below).  The 
solution  for  the  finest  mesh  of  135  x  58  grid 
points  differs  by  only  5%  in  stagnation  point 
location  from  the  92  x  58  mesh  solution.  The 
axial  distributions  of  other  variables  are  almost 
identical  between  these  two  meshes.  There¬ 
fore,  the  latter  mesh  was  adopted  for  most  of 
the  computations 


Findings 

Wake  Shock.  The  wake  shock  can  be  visu¬ 
alized  in  Fig.  5,  by  reversals  in  the  constant- 
density  contours.  Inspection  of  all  the  com¬ 
putations  shows  that  the  origin  of  the  shock 
lies  between  1  and  1.3  base  diameters  down¬ 
stream  of  the  base,  independently  of  Mach  and 
Reynolds  numbers.  The  location  of  this  point 
is  of  some  importance,  because  it  is  easily  rec¬ 
ognizable  in  photographs  of  free-flight  exper¬ 
iments,  and  it  might  be  used  as  a  reference 
for  other  features  of  the  wake  structure,  such 
as  the  location  of  the  rear  stagnation  point. 
The  wake  simulations  of  Reklis  and  Conti33 
at  Mach  6.32  indicated  no  correlation  between 
the  origin  of  the  wake  shock  and  the  location  of 
the  rear  stagnation  point,  because  the  former 
is  insensitive  to  Reynolds  number,  but  not  so 
the  latter.  The  present  investigation  confirms 
these  conclusions  for  Mach  4.32  and  6.32  but, 
as  indicated  below,  it  shows  that  at  Mach  10 
the  rear  stagnation  point  becomes  insensitive 
to  Reynolds  number,  and  its  axial  location  is 
indeed  close  to  the  origin  of  the  wake  shock. 

Rear  Stagnation  Point.  The  rear  stagna¬ 
tion  point  is  at  the  downstream  end  of  the 
recirculation  region,  and  it  indicates  the  ax¬ 
ial  extent  of  the  slow,  recirculating  flow.  The 
wake  flowfield  is  visualized  in  Fig  6  by  the  ve¬ 
locity  fields  and  streamlines,  for  selected  Mach 
and  Reynolds  numbers.  The  results  shown 
in  that  figure  for  Mach  4.32  and  6.32  suggest 
that  in  supersonic  and  low  hypersonic  flow  the 
position  of  the  rear  stagnation  point  is  sensi¬ 
tive  to  Mach  and  Reynolds  numbers:  increas¬ 
ing  the  Mach  number  moves  the  stagnation 
point  closer  to  the  base,  whereas  increasing  the 
Reynolds  number  moves  it  away  from  the  base. 
However,  at  Mach  10  the  location  of  the  stag¬ 
nation  point  is  insensitive  to  Reynolds  number. 
This  is  illustrated  in  Fig.  7,  which  compares  all 
the  present  results  with  the  measurements  by 
McLaughlin18  and  Blankson20,  and  the  corre¬ 
lation  by  Martellucci  et  al2i .  The  vertical  bars 
represent  the  uncertainty  in  the  McLaughlin 


and  Blankson  experiments  at  Reynolds  94.000 
and  86,000  respectively.  The  present  computa¬ 
tions  at  Reynolds  86,000  fall  within  these  un¬ 
certainty  limits,  and  are  considered  to  agree 
with  the  experimental  data.  The  disagree¬ 
ment  reported  by  Reklis  and  Conti33  between 
their  early  results  and  the  Blankson  data  has 
been  traced  by  the  present  investigators  to 
lack  of  sufficient  grid  resolution  in  the  early 
wake  simulations.  The  hypersonic  correlation 
by  Martellucci  et  al.  as  well  as  the  shape  of  the 
curve  faired  through  the  present  results  in  Fig. 
7,  suggest  that  at  Mach  10  and  beyond,  the  lo¬ 
cation  of  the  rear  stagnation  point  lies  at  ap¬ 
proximately  one  base  diameter  downstream  of 
the  base,  independently  of  Mach  and  Reynolds 
numbers. 

Flow  Separation.  Several  investigators 
( e.g .  Donaldson20,  Kayser27,28)  have  reported 
flow  separation  taking  place  beyond  sharp  or 
rounded  corners,  on  the  base  itself.  In  order  to 
investigate  this  phenomenon,  the  present  nu¬ 
merical  simulations  were  conducted  with  suf¬ 
ficient  grid  resolution  near  the  base  to  estab¬ 
lish  the  location  of  the  flow  separation  line. 
A  typical  base-flow  velocity  field  is  shown  on 
Fig.  8,  where  the  location  of  the  separation 
line  (indicated  by  an  arrow)  has  been  interpo¬ 
lated  between  two  mesh  points.  The  results 
of  all  the  present  computations  are  plotted  in 
Fig.  9  as  a  function  of  Mach  number.  There 
is  a  clear  trend  for  the  separation  line  to  move 
away  from  the  corner,  towards  the  axis,  with 
increasing  Mach  number.  At  Mach  4.32  and 
6.32,  increasing  the  Reynolds  number  moves 
the  separation  line  outwardly  towards  the  cor¬ 
ner.  At  Mach  10  the  separation  line  is  insensi¬ 
tive  to  Reynolds  number. 

Axial  Velocity.  Typical  plots  of  axial- 
velocity  distribution  along  the  axis  of  symme¬ 
try  are  shown  on  Fig.  10.  Near  the  base  the  ax¬ 
ial  velocity  goes  through  a  minimum,  beyond 
which  it  increases  asymptotically  to  its  far-field 
value.  The  location  of  the  minimum  shows  a 
weak  dependence  on  the  Mach  and  Reynolds 
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numbers,  and  it  correlates  well  with  the  axial 
location  of  the  center  of  rotation  of  the  recir¬ 
culating  How,  as  shown  graphically  on  Fig  8. 
This  correlation  is  plotted  in  Fig.  1 1 ,  where  all 
but  one  of  the  data  points  fall  within  a  :tl5% 
error  band.  The  good  correlation  is  attributed 
to  the  fact  that  in  the  low-speed  environment 
of  the  wake  the  rotating  flow  imposes  the  pres¬ 
sure  of  its  core,  and  the  corresponding  speed, 
on  the  wake  axis. 

Axial  Distribution  of  Static  Pressure. 
Typical  distributions  of  static  pressure  along 
the  wake  axis  are  shown  on  Fig.  12.  Two  fea¬ 
tures  of  the  pressure  distribution  are  the  min¬ 
imum  associated  with  the  velocity  minimum 
discussed  above  and,  in  some  cases,  an  over¬ 
shoot  of  the  static  pressure  beyond  the  free- 
stream  value,  which  it  should  approach  asymp¬ 
totically  in  the  far  field.  McLaughlin18  pointed 
out  that  the  overshoot  is  greater  at  lower  Mach 
number,  and  the  present  results  confirm  this 
statement. 

Base  Pressure.  Figure  13  shows  a  plot  of  the 
wall  pressure  at  the  center  of  the  base,  taken 
as  representative  of  the  base  pressure  Some 
radial  variation  of  base  pressure  was  observed 
in  the  numerical  simulations,  consisting  mostly 
of  a  decrease  in  pressure  from  the  center  of  the 
base  towards  the  corner.  The  computed  results 
show  considerable  scatter,  a  situation  that  is 
typical  of  other  studies  (logarithmic  scales  are 
widely  used  to  plot  base  pressure).  No  clear 
trend  with  Mach  or  Reynolds  numbers  emerges 
from  the  computations,  but  the  results  are  in 
general  agreement  with  the  extrapolated  data 
reported  by  McLaughlin  and  Blankson  ( loc . 
cit).  In  connection  with  extrapolating  to  the 
wall  the  pressures  measured  or  computed  in 
the  flowfield  to  obtain  base  pressure,  it  should 
be  noted  that  considerable  pressure  gradients 
exist  near  the  wall,  as  shown  in  Fig.  12. 

Conclusions 

•  Numerical  simulations  of  hypersonic 
wakes  show  good  qualitative  agreement 


with  experimental  data.  Major  flow 
features  show  good  quantitative  agree¬ 
ment. 

•  The  wake  flow  separates  on  the  base  of 
the  cone,  between  the  sharp  corner  and 
the  axis  of  symmetry. 

•  The  axial  velocity  has  a  minimum  near 
the  base.  Its  location  correlates  with 
the  location  of  the  core  of  the  recircu¬ 
lating  flow. 

•  The  static  pressure  on  the  axis  over¬ 
shoots  its  free-stream  value,  particu¬ 
larly  at  low  Mach  number. 

•  Mach  number  effects: 

Increasing  the  Mach  number  from  4  to 
10  has  a  significant  effect  on 
Shortening  the  recirculation  region 
Shifting  the  separation  line  towards 
the  axis 

And  a  weak  effect  on 

-Base  pressure 

-Origin  of  the  wake  shock. 

•  Reynolds  number  effects:  At  Mach  4 
and  6,  increasing  the  Reynolds  num¬ 
ber  lengthens  the  recirculation  region, 
and  shifts  the  separation  line  towards 
the  base  corner.  At  Mach  10,  Reynolds 
number  effects  are  negligible. 
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Fig.  3.  Topology  of  the  physical  domain  and  bound 
ary  conditions  for  wake  computations. 


Fig.  2. 


The  Reynolds  number  at  which  transition 
occurs  one  diameter  downstream  of  the  base, 
as  a  function  of  Mach  number.  Reference 
parameters  are  the  free  stream  and  the  base 
diameter. 


g.  6.  Velocity  field  and  streamlines  in  the  near  waKe,  for  selected  Mach  and  Reynolds  numbers.  The 
computed  location  of  the  rear  stagnation  point  is  labelled  “S";  experimental  values  are  indicated 
by  the  liars. 


Fig.  7.  Axial  location  of  the  rear  stagnation  point 
downstream  of  the  base,  normalized  by  the 
base  diameter.  Vertical  bars  indicate  the 
range  of  experimental  values.  The  horizon¬ 
tal  bar  represents  a  correlation  of  hyper¬ 
sonic  experimental  data. 
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Fig.  8.  Velocity  field  near  the  base,  for  Mach  6.32 
and  Reynolds  86,000.  Radial  (R)  and  axial 
(X)  coordinates  are  normalized  by  the  base 
diameter. 
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Fig.  9.  Radial  location  of  the  flow  separation  line 
on  the  base,  as  a  function  of  Mach  number, 
for  a  range  of  Reynolds  number. 
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10.  Effect  of  free-stream  Mach  number  on  the 
axial-velocity  distribution  along  the  axis  of 
symmetry. 


Fig.  11.  Correlation  between  the  axial  location  of 
the  core  of  the  recirculating  flow  (center  of 
rotation)  and  that  of  the  minimum  ar<al  ve¬ 
locity. 
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Fig.  12.  Effect  of  free-stream  Mach  number  on  the 
static  pressure  distribution  along  the  axis 
of  symmetry. 


Fig.  13.  Base  pressure  normalized  by  free-stream 
static  pressure.  Filled-in  symbols  indicate 
experimental  data. 
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